THE ENERGY LEVEL SHIFTS, WAVE FUNCTIONS 
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^ I We discuss the equations for the bound one-active elec- 

tron states based on the analytic solutions of the Schrodinger 
and Pauli equations for a uniform magnetic held and a single 
P^. attractive 6{r)-potential. It is vary important that ground 

p ^! electron states in the magnetic held differ essentially from the 

^ I analogous state of spin-0 particles, whose binding energy was 

I— !■ intensively studied more than forty years ago. We show that 

I ■ binding energy equations for spin-1/2 particles can be ob- 

>• . tained without using the language of boundary conditions in 

the 6-potential model developed in pioneering works. We use 
Q\ [ the obtained equations to calculate the energy level displace- 

^ ■ ments analytically and demonstrate nonlinear dependencies 

^ ■ on Geld intensity. We show that the magnetic held indeed 

^ . plays a stabilizing role in considered systems in a case of the 

^ ! weak intensity, but the opposite occurs in the case of strong in- 

>• I tensity. These properties may be important for real quantum 

^ ■ mechanical fermionic systems in two and three dimensions. 

^ ' We also analyze the exact solution of the Pauli equation for 

- - - an electron moving in the potential held determined by the 

three-dimensional 6-well in the presence of a strong magnetic 
held. We obtain asymptotic expressions for this solution for 
different values of the problem parameters. In addition, we 
consider electron probability currents and their dependence 
on the magnetic held. We show that including the spin in 
the framework of the nonrelativistic approach allows correctly 
taking the effect of the magnetic held on the electric current 
into account. The obtained dependencies of the current dis- 
tribution, which is an experimentally observable quantity, can 
be manifested directly in scattering processes, for example. 
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1 Formulation of the problem 

The effect of an external electromagnetic field on nonrelativistic charged particles systems 
(such as atoms, ions, and atomic nuclei) has long been investigated systematically (see, 
e.g., [1]-[I5] ). Although this problem has a long history, a set of questions still requires 
additional study. For example, a systematic analysis of the bound states of spin-1/2 parti- 
cles in an intense magnetic field has been lacking until now. We note that the basic results 
in the case of spinless particles were obtained using analytic solutions in nonperturbative 
mathematical treatments. As usual, the exact solutions of Schrodinger equations with 
Hamiltonians taking a particle bound by short-range potential in the presence of exter- 
nal fields into account are used. Furthermore, there is a rather common opinion that 
the role of the magnetic field in decays of quasistationary states is invariably stabiliz- 
ing [3], [H], [I5]. This view arises because the spinor states of electrons in an external 
electromagnetic field are usually neglected in nonrelativistic treatments, which is often 
inadequate [16], [17]. In this paper, we treat an essential part of these problems. 

We consider charged spin-0 and spin-1/2 particles bounded by a short-range potential 
{S potential) and located in an external stationary magnetic field with an arbitrary inten- 
sity. We note that the zero-radius potential is a widespread approximation for a multielec- 
tron atom field and especially for a negative ion field [12], [18]. Energy level displacements 
can be seen for the particle in a 5 potential and a magnetic field. The binding-energy 
equation is most appropriate for the investigating such states(|H], [12], [T9]). 

When an electron moves in a uniform magnetic field oriented in the z direction, the 
quantum mechanical system is invariant with respect to the z axis. The system then be- 
comes essentially two-dimensional in the xy plane. Many physical phenomena in axially 
symmetric quantum systems of electrically charged fermions (the quantum Hall effect [20j . 
high-temperature superconductivity [21], various film defects [22], etc.) can be effec- 
tively studied using the nonrelativistic equations of motion in 2+1 dimensions. A number 
of effects in constant magnetic fields, including certain types of doped two-dimensional 
semimetals, can be described using the Dirac equation in 2+1 dimensions [23], [23]. But 
there are many physical phenomena that occur in three-dimensional space [25], [26]. In 
this paper, we investigate the effect of a stationary uniform magnetic field on localized 
electron states in 2 + 1 and 3+1 dimensions (see also [27], [28]). 

The effects of external electromagnetic fields on bound nonrelativistic charged particles 
have already been studied systematically in detail over several decades (see, e.g., [T5] . 
[16], [29]). But although this problem has a very long history, several problems still 
need additional study. In particular, such problems include analyzing the effect of a 
strong magnetic field on bound charged particles with their spin states taken correctly 
into account. 
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There is a widespread opinion that consistently including the charged-particle spin is 
required only in studying relativistic effects. But this is not always the case for a strong 
magnetic field, as shown in [16], [30]. We note that efforts to analyze spin effects in the 
nonrelativistic approximation were made many times. Models precisely taking the effect 
of the field on the particle into account were used for this purpose. We emphasize that 
the wide experience in describing spinless particles bound by the (5-like potential in strong 
electromagnetic fields has been acquired precisely because the nonperturbative mathemat- 
ical approach was used. Developed methods based on the exact analytic solutions of the 
Schrodinger equation [TJ], [2S], [31] were also used to study processes taking particle spins 
into account. But because mathematical estimates are complicated in the case of spinning 
particles, the obtained conclusions are not always sufficiently convincing. In our opinion, 
the papers [27], [30], [32], where the effects of the field on the behavior of spinning and 
scalar particles were compared, are most correct from the standpoint of the necessity to 
take particle spins into account. For these purposes, the solutions of the Pauli equations, 
together with the solutions of the Schrodinger equation, were also used in the indicated 
papers. We here consider the effects produced by the action of a strong constant uniform 
magnetic field on an electron bound by a short-range potential. 

Under the action of the field, electron energy levels are shifted by a quantity determined 
from the transcendental equation for the energy. We note that this problem was analyzed 
in the case of a scalar particle in [31]. A similar problem for the electron in the magnetic 
field with its spin states taken into account was correctly solved in [27] quite recently. 

Our main purpose is to derive equations for the binding energy of a fermion in a field 
containing an attractive singular potential and a stationary external magnetic field in the 
two- and three-dimensional cases. We apply standard quantum mechanical methods using 
the expansion of the unknown wave function in a series at the eigenf unctions obtained for 
the fermionic system in the pure magnetic field. This method was used to study some 
physical examples of the effect of a constant magnetic field on charged particles bound 
by a single attractive 6 potential [33], [31]. This formalism differs in principle from the 
traditional derivation of wave functions in similar problems using the boundary condition 
typical for the 5 potential [1], [5], [H], [15] . 

It is very important that our approach permits developing a consistent investigation of 
the spin effects arising in an external magnetic field. An exact analytic expression for the 
wave function of a charged scalar particle in a state bound by the 5-potential and moving 
in a strong magnetic field was found in [31] using the Green's function of the corresponding 
Schrodinger equation. We note that the Green's function of a scalar particle in an external 
magnetic field was presented in the classic monograph [35] (also see [E]). As is known, 
the electron spin can be taken into account, for example, by passing to the nonrelativistic 
limit in the solutions of the Dirac equations describing the motion of a spinning particle 
in a given external field [36], [37] . 

The general structure of the paper and its main results are as follows. In the next 
section, we construct the equation for scalar particles with a low binding energy in a sta- 
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tionary external magnetic field based on an explicit solution of the Schrodinger equation. 
In Sec. Ill, we obtain the expressions for the energy of electron bound states in the S 
potential and an external magnetic field based on analogous analysis of explicit solutions 
of the Pauli equation. In Sec. IV, we discuss the equations for the binding energy of spin-0 
and spin-1/2 particles in the presence of the both weak and strong magnetic fields because 
particle spin was previously taken into account inadequately in similar problems. 

Furthermore, in Sec. V, we analyze the exact solution of the Pauli equation for an 
electron moving in the potential well determined by the three-dimensional (5-function in 
the presence of a strong magnetic field and to obtain asymptotic expressions for this 
solution in the case of different values of the problem parameters. In addition, in Sec. VI 
we consider the probability currents for the given particle and their dependence on the 
spin and magnetic field. 

2 A scalar particle in an attractive potential in the 
presence of a uniform magnetic field 

We consider a charge in a uniform magnetic field B specified as 

B = (0, 0, 5) = V X A, A = {-yB, 0, 0). (1) 
The Schrodinger equation in field ([1]) has the form 



d 

ih—il;{t,r)=ntp{t,r), v={x,y,z), (2) 
with the Hamiltonian "H is 

n- — (-ih— + —yX (3) 
2m \ dx c / 2m dy"^ 2m dz^ ' 

where m and e are the particle mass and charge. The particle wave function in field ([T]) 
has the form [19] 

where 

is the electron energy spectrum, u = \eB\/mc, and and Pz are the electron momenta 
in the x and z directions. 
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The functions 



are expressed in terms of the Hermite polynomials Hn{z), the integer n = 0,1,2,... 
indicates the Landau level number, = \/hc/\eB\ = ^yh/moo is the so-called magnetic 
length (see, for example, [25]) and yo = —cp/eB. 

We now study a simple solvable model. We consider the motion of a scalar particle 
in the three-dimensional case in a single attractive 5(r)-potential, where 6{t) is the Dirac 
delta function, in the presence of a uniform magnetic field. In fact, we must solve the 
Schrodinger equation 



1 
2m 



-ih— + —y] -f?—-h^—-h^5{v) 
ox c J oy^ oz'^ 



^E'{r) = E'^E'ir). (6) 



We can take solutions of Eq.(l6]) in the form 



oo „ 

*E'(r) = '^E'np^P^i^np^pAr) = Y^ / dp^dp^C E'np^p^^np^pA^) , (7) 



n,px,Pz 1=0 



where V'np^pz(i^) is the spatial part of the wave functions (jl]). 

The coefficients CE'np^p, can be easily calculated, and we then obtain the equation 



oo „ 
n=0 

where is a normalized coefficient independent of the field and 

. I E p^^ 

A = h — — . (9) 

2 Hlo 2mhjjj 

Integrating over p^ gives ([8]) in the form 

oo ^ 



n=0 



-f A)V2- 



It is easy to see that Eq.f llOp implicitly defines the energy of a bound localized electron 
state in the magnetic field. We note that (11 01) is consistent with analogous result in |34j . 
where this equation was solved numerically. But Eq. ( ITO|) can be analytically reduced to 
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a simpler form. Indeed, we can sum over n in the right-hand side of Eq.f llOl) using the 
representation 

oo 

' dt. ill) 



[n + A + iefl^ J 



As a result, Eq. ffTOl) becomes 



e 4 / e '"^ 



1 = iViV/ia; — ^ / —j:: r^dt, (12) 



where A'^i is a real constant independent of the field. Because the required energy 

E' = -\E'\, (13) 

must be negative, we can rotate the integration contour through the angle 7r/2 in the 
complex t plane. We thus obtain the real expression 

oo ^ 

2^/^ J tV2sinh(t/2) ' ^ ' 



where E' = > 0. If we eliminate the magnetic field, then ( HM takes the form 

oo 

r p-Eot/h 

- J -^dt^ (15) 



where Eq = \Eq\ is the absolute value of the binding energy of the particle in the 6- 
potential without the action of the external field. Subtracting f|T5l) from f fT^ and removing 
the integral divergences in the lower limit by the standard regularization procedure, we 
obtain 

oo oo 

Eot/h _ -Et/h p -Et/h / \ 

dt= -^-Adt, (16) 



t3/2 J t3/2 ^sinh(ait) 





where ai = ^. From (1161) . it is easy to obtain 



oo 



2y/7r J X'^i^ \smh(ax) 
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where a = f^, which is consistent with the analogous equation obtained by the well-known 
method using boundary conditions of wave functions in the 5-potential model [5].|14j.[32]. 
Expanding the integrand function in ( fT71) in the weak- field limit hu <^ 2£'o, we obtain 

V 48 Eo^ 576 Eo^ 

We note that the quadratic term in (fT8|) coincides with analogous result in [5]. 

To consider the strong-field case hj > 2Eq, we reduce the right-hand side of (fT7|) to 
the analytic form 

1 1 /I 1 1 \ 

:C o.o + 7T-h (19) 



^ 72^2' 2 2a 



where oq = and Cl'^iP] is the Hurwitz (generalized Riemann) zeta function. The 
vahdity range of (1191) can be found somewhat wider that was initially assumed. In deriving 
(lTi|) we assume that E' < 0, but we can see from (fT9|) that argument of zeta function can 
continuously reach the values 

1/2 + l/2a > 0. 

This condition limits the required binding-energy spectrum by 

E' < — . (20) 

The physical meaning of this condition is the restriction to the continuous spectrum of the 
scalar particle in the magnetic field by the value of fl20l) (i.e., E' > tvjj/2). We note that 
after a change of variables in ( fT9l) . it is coincides with the basic equation in [5], where the 
case of scalar particles in the magnetic field was considered and an analogous conclusion 
about limitation of continuous spectrum was drawn. 
Expanding Ci^iV) for p ^ 1 gives 

C(1/2,P) = ^ + C(l/2) - ic(3/2)p + ^C(5/2)/ + Ob]l (21) 

Substituting fl2T|) in f fT9|) . we obtain the explicit equation for the bound-state energy in 
the strong-field limit 



E' = hw\ 0.205 - 0.452W ^ - 0.367 ^ ) . (22) 
\ \ rkj rkj I 

We emphasize that in superstrong magnetic fields, expansion (l22l) gives the upper limit 

E' = 0.205 
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for the binding energy of the scalar particle, which does not contradict condition (!20|) . 
Furthermore, it can be seen that this limiting value is independent of the particle energy 
in the absence of the field and is completely determined by the magnetic field intensity. 

It is interesting to compare the obtained results with the results in the case of two- 
dimensional model. The analogue of Eq. lfTOl) in the two-dimensional case is 



oo 



]^ r -Et/hw 

1 = — / ^^dt, (23) 

Stt J sinh(t/2) ' ^ ' 



which coincides with the corresponding result in [33]. But our regularization procedure 
here essentially differs from [33]. As before (see ( 1161) ). we remove the magnetic field and 
obtain 



oo 



1 r p-Eot 


With a simple calculation similar to that in the three-dimensional case, we can write 

oo 

ln#= [—(-^-^)dx, (25) 



Eq J X \sinh(aa; 



where a = hw/{2E) as before. In the weak-field limit, we obtain 



E = E,{l- ] (26) 



2AEc 
from (|25|). 

To consider the range hw > 2Eq, we must first calculate the integral in the right-hand 
side of Eq. (l25l) analytically. 



- (I) - '"P") ^ * (^) ■ '''' 

where "^{x) is a logarithmic derivative of Euler gamma function. We then have the basic 
equation in the two-dimensional model 

-ln(2.„) = *(i + i-), (28) 

where Oq = huj/{2EQ). 
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In the strong-field limit, after evaluation the function 



^(p) = - C + + ^vl/(2)(iy + + o[p]^ (29) 
p b 2 90 

where 

Jn 

¥-\z) = —^{z) = {-ir^\K{n + l,z), 

we can write f l28|) as 

where C = 0.577... is the Euler constant. 

The solution of Eq. (!30|) . which explicitly determines the bound-state energy, can be 
written as 



E' _ 1 6(C - \n{hu/Eo)) + ^I^tx'^ + 36(C - ln(na;/Eo))^ 
^ ~ 2 2^^2 

For ln(^ci;/-Eo) ^ 1 we obtain 

E' \ 1 C (t"'^ 



(31) 



+ . 3,, ; '^ + 0[ln(na;/£;o)]^ (32) 



/kJ 2 ln(na;/^o) ln^(^/^o) ln'^(/iu;/£'( 

from ( l3Ti) . Considering the properties of ^{^z\ we see that expansion ( l32l) is correct for the 
binding energy E' =< fiLo/2. Furthermore, this limiting value, as before, is independent 
of the particle energy in the absence of the field. But there is an essential difference from 
the three-dimensional case. For supperstrong magnetic fields (when not only hw/Eo is 
large but also ln(/za;/£'o) >> 1), the upper limit of the shifted binding-energy level in the 
considered model tends directly to the boundary of the continuous spectrum. 

3 A spin particle in an attractive potential in the 
presence of a uniform magnetic field 

It is vary important that we can use the present approach to study the spin effects in 
magnetic fields in the same way. The case of a spin-1/2 particle can be calculated based 
on exact solutions of the Pauli equation. The Pauli equation in the field ([1]) has the form 

d 

ih—ip{t,r)=ntp{t,r), r={x,y,z), (33) 
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with the Hamiltonian "H is 

H- — i-ih— + — - — — - — — + ^^a B (34) 
2m \ dx c J 2m dy"^ 2m dz"^ ' 

where /i = \e\fi/2mc is the Bohr magneton, m is the mass of a electron and 



1 
-1 



is the 2;-component of Pauh matrices. The last term in fl34j) describes the interaction of 
the electron spin magnetic moment with the magnetic field. The electron wave function 
in field ([T]) has the form 

^np,p,.(i, r) = ^V'np,p,(t, r) J ]^ J ^ , (35) 
where 'ipnp^pX'^^'^) is the solution of the Schrodinger equation in field ([T]) (see(Ilj)), 

Ens ^ (^'^ + + 2m sfuJ^ (36) 

is the electron energy spectrum, uj = \eB\/mc, s = ±1 is the conserved spin quantum 
number, and Px and Pz are the electron momenta in the x and z directions. 

It is very important that the electron ground state in a magnetic field differs essen- 
tially from the analogous state of spin-0 particles. Moreover, the continuous spectrum 
boundaries differ for spinor and scalar particles. For example, if the continuous spectrum 
of a scalar particle begins from E' > hw/2, then the continuous spectrum of an electron 
begins from E' > hoj for the spin directed along the magnetic field direction and from 
E' > for the spin directed against the magnetic field direction. 

Taking the interaction of the electron spin magnetic moment with the magnetic field 
into account, we can write the energy equation in the three-dimensional case in the form 



E 



'En 



2v/? J t3/2 \^sinh(at) 



ate 



sat 



1 dt. 



(37) 







were s = ±1 respectively corresponds to spin orientations along or against the magnetic 
field direction. Expanding the integral in f l371) for a << 1, we obtain the equation 



'En 



shu ^/e ( hhj 



4^ 12 V2^( 



(3J 
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The solution of Eq. (!38|) has the form 




E 

Eq 



2 ( -12Eo - V^^/A8Eo'^ + sih^u^Eo - AShuEo^ 



(39) 



Expansion ( l39i) can be written in the weak-field limit as 



E 



1 - s 



hjj 1 h^u}"^ 



2En 48 E^ 



(40) 



It is easily seen from ( HOl) that the energy level Eq = —\Eq\ existing in the 5-potential 
without a perturbation is shifted under the magnetic field action by +^ (upward) in the 
case s = 1 and by (downward) in the case s = —1. But the depth of the arrangement 
of energy levels with respect to the continuous spectrum boundaries is the same in these 
two cases. We note that it has the same depth in the case of spin-0 particles. 

Integrating of the right-hand side Eq. (!37|) . we obtain the equation in the analytical 
form 



1 



1 1 
2' 2 



s 
2 



1 

2a 



(41) 



In the strong-field limit hw > Eq, we can write the Hurwitz zeta-function in Eq.f l41l) as 

1 1 , c[i/2] c[3/2] n 



1 1 
2' 2 



1 
2a 



V2 



l+s I E_ 
2 ^ hid 



Finally, we write the Eq.( l4Tl) as 



V2 2V2 




'0 I C[l/2]^ C[3/2]^3 

hi^ V2 J 2V2 



(43) 



where 



X 



l + s E 
2 ^ hw' 



The solutions of Eq. (H3!) for different spin values s = 0, +1, —1 can be represented as 



(s I E E \ 

0.205 + -- 0.452W-^- 0.367-^ . (44) 

In Fig. 1, we show graphs of solutions of Eq. fBTl) for different particle spin values with 
£^0 = 1 and hu = IOO-Eq, which are typical values for experiments on the bound-state 
energy in semiconductors at low temperatures in magnetic fields H ~ lO^G. 
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Figure 1: Solutions of Eg. (14 ip for different particle spin values. The horizontal straight 
line is y = ao~^^^. The curves are y = — 2^^/^C[l/2, (1 + s)/2 + l/2a] and from left to right 
s = —1, 0, 1. The parameter W = E' / Eq is also shown. 

It is obvious that the approximate solutions of dH]) are quite near the intersections of 
the graphs of left-hand side of Eq. fHTl) (the straight line in Fig. 1) and of the right-hand 
sides of Eq. (l4ip (the curves in Fig. 1). We emphasize that the dependence of energy level 
shifts on the particle spin does not disappear in the strong-field limit. Moreover, the 
continuous spectrum boundaries are shifted in the cases s = and s = 1. Hence, the 
perturbative displacements of the binding-energy levels (as in the weak-field limit) are at 
the same distances from the continuous spectrum boundaries in all cases. 

We now consider spin interactions in the two-dimensional case. According to our 
approach, we can write 



where the particle spin direction, as before, is represented by s = ±1. In the weak-field 
limit, we obtain 



from Eq. (H5|) . To consider the strong-field limit, we must calculate the integral in Eq. fl45p 




(45) 








(46) 
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for s = ±1 in analytic form 



oo 



e Uaxe _ A ^ -2a^^^ - ln(2a) - ^ ( - ) • (47) 



X \sinh(ax) / 2 \2a 



We can then write the equations for energy displacements for hu > Eq as 

For the case s = — 1 in the strong- field limit (In ^ >> 1) we immediately obtain 

E - -hu I ^ + J 1 (49) 
from Eq.(jl5]). 

For the opposite spin orientation (s = 1) in Eq. (l47l) . we must first use the recurrent 
relation for \l/(p), 



1 + ^(x) = $(1 +x). (50) 

X 

Using the asymptotic expansion for \l/(p) (see (1291) ). then obtain 

E' = huj 11- - 1 . (51) 

The dependence on the spin parameters can be interpreted as in three-dimensional 
model. We can write this equation in the form 

,„f^^._,,fi±l+ IV (52) 



EoJ V 2 2a 



In Fig.2, we show graphs of the solutions of Eq. fl52l) for different particle spin values 
with Eq = 10~^eV and hw = IOOEq. The main difference from the three-dimensional 
model is that the perturbative binding-energy levels converge to the continuous spectrum 
boundaries in a superstrong magnetic field in this case. 
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Figure 2: Solutions of Eg. (15 2 p for different particle spin values. The horizontal straight 
line is y = Infuj/Eo. The curves are y = — ^E'[(l + s)/2 + l/2a] and from left to right 
s = —1, 0, 1. The parameter W = E' / Eq also shown. 

4 The binding energy of spin-0 and spin-1/2 parti- 
cles in the presence of the both weak and strong 
magnetic fields 

We have shown that the effect of a magnetic field on localized electron states leads to 
equations for the binding energy of spin-0 and spin-1/2 particles. In the weak- field limit 
{tvjj <^ Eq), the energy displacements of scalar and spinor particles are described by the 
expressions 



s = 



s = 



s = 



= 
= 1 
-1 



2Eo 
nw 

Eo 



+ 
+ 




(53) 



in the three-dimensional case 



and 



s = 



s = 



s = 



= 
= 1 
-1 




(54) 



in the two-dimensional case. From Eq.( l53l) and Eq. (l54p we can also obtain 
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= 0: E = Eo-^, 

1: E = Eo-- 



2 245^0 

, = £ = + (55) 

where 5 = 2 and 5 = 1 in the three-dimensional and in the two-dimensional cases corre- 
spondingly. 

The dependence on the particle spin does not disappear in the strong-field limit 
{J\bj ^ E'o). In three-dimensional case, the perturbative energy levels approach specific 
spectral values determined by the magnetic field intensity. For different particle spin val- 
ues, the displacements of the binding-energy levels as before [see Eq. fl53|) and Eq. flMl) ] are 
at identical distances from the continuous spectrum boundaries and can be represented 
in the form 



s = 


2 


- E' 


s = 1 




- E' 


s = -1 




- E' 


Removing the braces 


in Eq4 




s = 


: E' 




s = 1 


: E' 


s = 


-1 : 


E' = 




0.367^ 1 . (56) 



E' = 0.705hu - 0A52^/E^ - CSGTE-o 



(}.295huj -OA52^/EohuJ -0.3Q7Eo. (57) 

In particular, from Eq.( l57p one can immediately see that the value of the binding 
energy level is positive both in the case of a spinless particles (s = 0) and in the case 
where the electron spin in oriented along the magnetic field direction (s = 1), but it 
remains negative in the case where the electron spin is oriented against the magnetic field 
direction (s = —1). 

It can be easily seen that the dependence on the spin parameters in the two-dimensional 
case can be written analogously, 

1=1 I - E'\=^ + -^, (58) 
. = -1 - E' 
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and without braces we have 




C hw 



s = 



: 




) 




) 



1 : 



E' = hu- 



c hw 



s = 




) 




) 



E' 



c hu 



(59) 



s = 



-1 : 




) 




) 



It can be seen from Eg. fl58p and Eq.f l59p that the energy levels in the basic terms 
are also independent on the particle energy in the absence of the magnetic field. The 
distinctive feature of this case is that the binding-energy levels for superstrong magnetic 
fields [when ln(^/£^o) >> 1] directly approach the continuous spectrum boundaries for 
all considered spin values. 

We have shown that the energy levels of a polarized electron under the action of a 
weak magnetic field for different particle spin values are shifted similarly in the three- 
dimensional and two-dimensional models. We have the line displacements as the levels 
themselves for s = 1 and s = — 1 and analogous shifts of the continuous spectrum bound- 
aries of for s = 1. We also have the same picture in the case of a spinless particle with 
the line shift of the continuous spectrum boundary. Clearly, in case of weak intensity, 
a magnetic field indeed plays a stabilizing role in the considered systems because the 
depth of the perturbative binding-energy levels from the continuous spectrum boundaries 
are shifted downward under the field action independently of the particle spin. But our 
results show a nonlinear dependence on the field intensity in the strong-field limit. Nev- 
ertheless, the continuous spectrum boundaries in the cases s = and s = 1, as before, 
have a linear dependence on the field in this limit. In superstrong magnetic fields, the 
binding-energy levels can approach the continuous spectrum boundaries. The distinctions 
can be formulated as follows. In the three-dimensional model, there is a fixed depth of 
the energy levels from the continuous spectrum boundaries that is the same for all spin 
values. In the two-dimensional model, the energy levels in a superstrong magnetic field 
tend asymptotically to the continuous spectrum boundaries. But in both cases, the sys- 
tem instability increases in strong magnetic fields. This conclusion therefore disproves the 
opinion that a magnetic field always plays a stabilizing role in systems of bound particles. 
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5 The exact analytic solution of the Pauli equation 
for the attractive three-dimensional (5-well and its 
asymptotic expressions 

The Green's function obtained in [36] is the solution of the Pauh equation with a 5-source 
and can be represented as an integral over time [37], [28]. This integral determining the 
Green's function admits a Wick rotation to the lower complex half-plane of the variable t 
(see, e.g., [25]). This operation makes the integral purely real. As a result, the stationary 
Green's function can be represented in the form 

Gw(r,0)= ^e—-»^y/'x 

where m is the particle mass, = eH/mc is the cyclotron frequency, e is the absolute 
value of the particle charge, H is the strength of the uniform magnetic field oriented 
along the z axis, s = ±1/2 is the spin number, and the argument of the exponential in 
the integrand is in fact the classical action function 

S = - \mu:y coth(^) + {W - s ■ hwjt (61) 

(although the argument of the exponential formally contains the Planck constant, this 
dependence vanishes because of the shift in the energy of a bound spinning particle in a 
magnetic field [27]). In formulas (1) and (2), x, y, and z are the Cartesian coordinates, 
= x'^ + y'^, and W < is the total energy of the bound particle. We note that (1), (2) is 
the ordinary propagator of a charged particle moving in a magnetic field and is continued 
to the range of negative energies. 

We write the spatial part of the wave function in the form 

^(f ) = Ne-'^'^H-y/^ A sinh-i e^/^ (62) 

where is the normalizing coefficient. We pass to more natural variables using character- 
istic scales of the problem. As an energy scale, we take \ Wq\, which is the absolute value of 
the particle binding energy in the absence of a magnetic field. We let w = W/\Wo\ denote 
the dimensionless binding energy in such units. It is equal to —1 in the absence of the 
external field. Of course, w depends on both the external field H and the spin direction in 
the general case [27] • It is convenient to measure the magnetic field in the dimensionless 
units h = huj^/\Wo\. The de Broglie wavelength Iq = h{2m\Wo\)~^^'^ of the particle for the 
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zero magnetic field can serve as a natural coordinate scale, and the quantity to = ^|W^o|~^ 
can serve as a time scale. Using such units for the spatial part of the wave function, we 
obtain 

= r jfrexp[-SV(4.)-;,^^coth(/W2)/8- ar]_ ^^^^ 

Jo T ' 1 — e 

where N is the normalizing coefficient in the new variables and x, y, z, and r are the 
dimensionless coordinates and time. The analogy between the integrand in (4) and the 
classical Planck formula for the black-body emission spectrum is interesting. The quan- 
tizing character of the magnetic field (as the quantum character of radiation) is not man- 
ifested until the exponent factor hr in the denominator in the right-hand side of (4) is 
sufficiently small. We also call attention to the phase factor g-^'^^S/^^ whose presence 
explicitly demonstrates the existence of the orbital probability current. 

To calculate the normalizing coefficient in (4), we must evaluate several integrals. The 
integrals over the coordinates are assumed to be purely Gaussian, and the integral over 
time can be reduced to the generalized Riemann zeta function C(3/2, e) by a simple change 
of variables. As a result, we obtain 



2nlf" \2' h 

where w = Wsh—w is the absolute value of the so-called effective energy of the particle. The 
parameter w appears because the lower edge of the continuum, together with the bound 
states w, is shifted in the magnetic field. Consequently, the binding effective energies are 
now measured from new boundaries determined by the expression Wsh = {s + l/2)h |27j . 
It also follows directly from this relation that the shift Wgh in the continuum depends 
explicitly on both the magnetic field and the particle spin orientation. But it was shown 
in [27] using weak {h <^ 1) and strong {h ^ 1) fields as examples that w is independent 
of the spin. 

It is easy to see that the integral in (4) can be related to Laplace- type integrals [38], |39j . 
The relevant integration domain is determined by a neighborhood of a single point of the 
exponential maximum. We consider the contribution of the saddle point to the integral 
in (4) written in the form 



g[u) = h — coth I 



/o Mi/2i_e-n' ^ 4m 8 \2 J 

Here, u = hr, and the saddle point is a root of the equation 



- + I^^- (65) 



Bu^-A = T-—-, A = Az'^p-^, B = IQwr'^h"^. (66) 

sinh (m/2) 
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Figure 3: The solution of transcendental equation (7) in the plane z = 1.6 for p = 1. The 
solid curves correspond to the left-hand side of Eq. (7) for different values of the magnetic 
field, and the dot-dashed curve is the graph of the right-hand side of Eq. (7). 



A graphic illustration of the search for the solution of this transcendental equation is 
shown in Fig. 3, which shows the family of parabolas corresponding to the left-hand side 
of Eq. (7) for different values of the parameter h (the solid curves) and the graph of the 
function in the right-hand side of the equation (the dot-dashed curve). It is easy to verify 
that the parabola branch can intersect the monotonic function sinh~^(u/2) in various 
ranges of the integration parameter u. Hence, in the range u <^ 1, the right-hand side of 
the equation differs slightly from the constant, and the saddle point, which is the root of 
the quadratic equation in this case, is given by 



Uq 



2Vw 



(67) 



In the other limit case u ^ 1, the right-hand side of Eq. (7) almost vanishes, and for the 
saddle point, we have 

hz , , 

uo = (68) 



Finally, for the intermediate range, the solution of Eq. (7) can be written in the approxi- 
mate form 



Mo 



6 + 6 



5W 



1 + XT^+ ^ 



~Z^h^ 



24w Vlw 2304M72 



1/2. 



1/2 



(69) 
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It is easy to see that this solution in the hmits of weak {h <^ 1) and strong (/i ^ 1) 
fields respectively transforms into (8) and (9). In this case, the evaluation of integral (6) 
obtained by the saddle point approximation is written as 

g(«o/2)-/it/(uo)_ (^YO) 



We now turn to studying the wave function in the limits of weak and strong fields. 
We see in what follows that these concepts require some more accurate definitions in the 
problem under consideration. If ^ 1 and if p and z are not very large, then solution (8) 
reduces to 

uo = —j^ <C 1, (71) 



where f = \/z^^p^ is the dimensionless radial coordinate. In the limit under consider- 
ation, this expression in fact determines the spherical type of symmetry of the electron 
cloud. 

As already stated in the introduction, the behavior of the bound energy level of a 
spinning particle located in a (5- well in the magnetic field was studied in [27]. It is impor- 
tant that the obtained solution has stable asymptotic expressions with respect to the spin 
variable in the limits of weak and strong fields, i.e., the effective energy characteristic w is 
independent of the spin in both cases. Although the computational technique used in [27] 
differs slightly from techniques used for scalar particles in the previous papers [H], [31], 
the results obtained in the limits of weak and strong fields in [27] agree completely. In 
particular, according to [27], we have u? ~ 1 -|- /i/2 for small magnetic fields, whence it can 
be directly seen that there is no dependence of w on the particle spin. 

Taking the foregoing into account, we use the estimate for the zeta function 

which is applicable in the range h <\. In particular, the graphs shown in Fig. 4 demon- 
strate that such an asymptotic expression is valid. Therefore, for the range under consid- 
eration h < \ (and not very far from the 5- well), we obtain the expression 

2 /tF ( % \ 
ip{f)^N ^_ expl—-hxy — fj (73) 

for the spatial part of the wave function. 

To represent the results graphically, it is convenient to use the function x(^) = rtp{f) 
having no singularities at zero and determining the spatial distribution of the probability 
density of the electron cloud in the spherical system of coordinates: 

dW ~ 27r|x(r)|2sin(^) dO dr. 
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Figure 4: The exact dependence of the function ({3/2, w/h) on the magnetic field h (the 
dashed curve) and the dependence determined by the right-hand side of formula (13) (the 
solid curve). 

Figure 5 shows the graph of |x(^)P- this case, the de Broglie wavelength is assumed to 
be unity. We stress that this estimate for the wave function was obtained in the vicinity 
of the (5-well in the weak field limit. But we note that the field of the order of several 
tenths of the interatomic field is not weak in the ordinary sense. 

Far from the 5-well, even in the weak field approximation, the case changes cardinally. 
For /i5 S> 1, from formula (10), we have 

uo = > 1 (74) 

for the saddle point. As a result, the asymptotic expression for the spatial part of the 
wave function becomes 

This expression is important for a qualitative analysis. The solution has an axial symmetry 
in this range. It is clear that the value of the field determines the range where the spherical 
symmetry typical of the bound s state transforms into the axial symmetry inherent in the 
wave functions of particles in a purely magnetic field. In addition, it can be seen that 
at large distances from the (5-well even in the weak-field case, only the ground level of 
the effective energy makes the main contribution, i.e., the level located at the smallest 
distance from the continuum edge [1 + h/2 in our adopted units). This agrees completely 
with the conclusions in [31], [27]. In the case where the processes with free particles in 
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Figure 5: The dependence of the dimensionless squared function on the transverse 

(x) and longitudinal (z) coordinates for fields h < 1. 



an external magnetic field are considered, it is usually assumed that the weak field does 
not always exhibit its quantizing character, yielding only small corrections to the cross 
sections of the corresponding processes. In contrast, it is traditionally assumed that only 
several minimum-energy levels contribute in the strong- field case [25] . 

In contrast to this, in the case of a bound state, any arbitrarily weak field behaves as a 
strong field at large distances from the center. In particular, it follows from (16) that the 
magnetic field compresses the electron cloud not only in the direction transverse to the 
field (which can be expected) but also along the field. But if the shift in the ground level 
of the effective energy in the magnetic field is neglected, i.e., if it is assumed that w = 1, 
then such effects cannot be described. We note that under these simplifying assumptions, 
solution (16) agrees with the basic function used in |33j . 

For a strong field, the case is more complicated. The reason is obvious because the 
intermediate range in which the symmetry transforms from the spherical into the axial one 
in this case belongs to the range where the wave function differs significantly from zero. 
Using the results in [27] on determining the effective energy of the electron in strong fields 
{h 3> 1), w ~ 0.295/;., we consider a small neighborhood of the (5-center. For f -C l/Vh, 
the position of the saddle point is determined by expression (12). For the wave function. 
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using relations (4) and (11), we obtain an analogue of expression (14) but with a different 
effective energy: 

fV0.295/i 



(76) 



Therefore, supplementing our conclusions on the behavior of the function at large distances 
in a weak field, we note that any arbitrarily strong field behaves as a weak one in the 
direct vicinity of the attractive (5-center. The reason is that the depth of the (5-well is 
much greater than any finite shift in the energy level in the magnetic field. The equation 
for the bound level energy w{h) follows precisely from the fact that the character of this 
asymptotic expression is independent of the external magnetic field [H], [31], |27j . 



K = 2000 




Figure 6: The dependence of |x(^)P ^ the axis p = for the field h = 100. The 
points correspond to the calculation using formula (4), and the solid curves (successively 
from the bottom up) determine the strong field approximations: the contribution of the 
ground level (18) and the contributions for A' = 20 and K = 2000 calculated using 
formula (20). 

The wave function becomes axially symmetric in the strong field at distances z ^ 
l/y/h. For h ^ 1 and fixed p ^ and z ^ 0, we obtain formula (15). For the wave 
function, from relations (4) and (11), we obtain an analogue of formula (16) in the case 
of strong fields: 



N 



TT 



0.295/1 



exp (^-^hxy - ^ - zV0.295h^ . 



(77) 



It is easy to see that also in this case, only the ground level of the effective energy 
contributes. 
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Spherically symmetric estimate (12) and axially symmetric estimate (15) of the max- 
imum point agree in a neighborhood of the straight line p = 0. It is interesting that the 
dependence of the bound level energy on the magnetic field was first obtained in |31j . 
where the limit transition in the expression for the wave function was used precisely along 
this line. Therefore, to extend the applicability range of the obtained expression (18) to 
the range of small 2, we can expand the denominator in initial formula (4) into a geometric 
series in terms of partial effective levels, 



replacing the Landau levels in a purely magnetic field [19], and use estimate (18) for each 
term. The wave function in this case becomes 



It is easy to verify that if the maximum number K of levels taken into account increases, 
then the applicability range for the obtained formula can be extended to arbitrarily small 
values of the coordinate z. 

In Fig. 6, the squared function |xp calculated using exact formula (4) is compared 
with asymptotic expressions (18) and (20) for the field h = 100 and the axis p = 0. The 
larger the parameter \/hz is, the better the series in solution (20) converges (and the first 
term of series (18) consequently approaches the exact solution). This turns out to be 
important because the probability current is maximum in this range. 

Figure 7 shows the comparison of the expressions for |xp in the case of three different 
values of 2 (2 = 0,0. 3, 0.6) and the strong field h = 100 calculated using exact formula (4) 
and approximate formula (11) with (10) taken into account. 

We briefly summarize the results in this section. The solution of the Pauli equation is 
spherically symmetric for a weak magnetic field in the entire relevant range of coordinates. 
It is in fact determined only by the character of the singularity in the (5-well; all effective 
partial levels Wk contribute to the formation of this singularity (see formula (19)). The 
quantizing character of the magnetic field (i.e., the tangibility of the contribution of indi- 
vidual energy levels) is manifested only far from the center in the region where the wave 
function is exponentially suppressed. For the strong field, the character of the solution 
symmetry changes in the relevant range of coordinates. The relative contribution of indi- 
vidual levels for a fixed field strength is mainly determined by the longitudinal coordinate 
(in the magnetic field direction), although the electron cloud is mainly compressed by the 
magnetic field in the transverse direction. 



Wk = w + kh, /c = 0, 1, 2, . . . , 



(78) 




^~zy/{0.295+k)h 



(79) 



y/{0.295 + k)h' 
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Figure 7: The dependence of |xP on p for the field h = 100 and different z. The points 
correspond to the calculation using exact formula (4), and the solid curves correspond to 
strong field approximations using formulas (10) and (11). 

6 Probability currents of a particle bound by the S- 
potential in a magnetic field 

The general expression for the nonrelativistic probability current of the spinor was given 
in the classic monograph [19]. Using the usual tensor notation, we represent this current 
Jfe in the form 

where = 1,2,3 are the spatial coordinate indices, a = 1,2 are the spin indices, Vfc is 
the usual gradient operator, Ak is the vector potential of the external magnetic field, /i is 
the electron magnetic moment, Ekpq is the unit totally antisymmetric tensor, and are 
the binary sigma matrices. The complete normalized solution of the Pauli equation can 
be represented in the form 

where the spatial part ilj{f) is given by (4). The first and second terms in (21) determine 
the gradient (orbital) current, and the last term corresponds to the spin current. 

Contracting spin indices with the explicit form of the sigma matrices taken into ac- 
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count, we obtain the spatial components of the total probability current 



2m V dx 



ijj — ijj' 



dx J mc 



e . _ , /ic s 9 , _ 

+ — AxV- V - — — — V'), 



e |s| dy 



2m \ ay ay/ mc e |s| ax 



2m \ az oz J mc 

We must fix the gauge of the vector potential for further calculations. Choosing it as 
= —Hy and Ay = = and passing to the dimensionless coordinates, we obtain 



J 



\Wo\ 
2m 



1/2 



h 



ip*ip{jx - ly) + 



{2m\Wo\f/^ lie s 



2h 



e \s\ 



ip*tp0x-ly), 



where i and j are the unit vectors of the Cartesian system and the function ip^f) has the 
form 



-hr 



dr 1 + e 



ri/2 [1 _ e-'^-]' 



exp 



h 



—p cothi — — wr 



At 8 



\2 J 



Assuming that the electron magnetic moment is equal to the Bohr magneton, we obtain 



J 



\Wo\ 
2m 



1/2 



^^*(^ + 2siIj){jx - iy). 



(82) 



Using exact expression (4) for the spatial part of the wave function, we obtain the 
formula for the probability current in the case of an arbitrary spin orientation and field 
strength: 

J = M/.n-^(|f)(j5-?y)x 

dr exp[-57(4r) - hp^ coth{hT/2)/8 - wt] 



oo ^ 

71/2 I _ g-hr 

°° dr exp[-5V(4r) - hp^ coth{hT/2)/8 - wt] 



T 



1/2 



1-e 



X 



l + 2s-coth( — j 



(83) 



where M = mlWol"^ /Att'^H^ is a dimension factor independent of the magnetic field. 
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We obtain the estimates for the current in the weak-field approximation. In this case, 
we have 

ip(r) ^ ■ —^^(1 + f) exp ( —-hxy — r] . 

Also taking (14) and (23) into account, we write the expression for the probability current 
in the weak magnetic field: 




Figure 8: The dependencies of ln(|J|/Af) on the transverse coordinate in the planes 
z = const for the weak field /i = 0.1 in the cases of (a) the particle with spin directed 
opposite the field and (b) the scalar particle. The points correspond to the calculation 
using exact formula (24), and the solid curves correspond to weak- field approximation (25). 

Figure 8a shows the dependencies of the logarithm of | J|/M on the transverse coordi- 
nate p in different planes z = const for the field h = 0.1; they are calculated using exact 
and approximate formulas (24) and (25) for s = —1/2, i.e., for the particle with the spin 
oriented opposite the field. In the case of the particle with the spin along the field, only 
the orientation of the vector J changes, i.e., the direction of the particle rotation changes. 
This occurs because the second term in (25) in the vicinity of the (5-well turns out to 
be much greater than the first term. In other words, the orbital probability current in 
this case is negligibly small compared with the spin one. They become comparable only 
at a distance r ^ 1/h from the center, where the current in the weak field {h -C 1) is 
exponentially suppressed. This conclusion is confirmed by comparing Fig. 8a with Fig. 8b, 
where the same quantity ln(|J|/M) is constructed for the scalar particle. In this case, 
the spatial part of the wave function in the weak field (see formula (14)) is independent 
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Figure 9: The dependencies of In (| J|/M) on the transverse coordinate calculated using 
exact formula (24) in the planes z = const for the strong field h = 100 in the cases of 
(a) the particle with the spin directed opposite the field and (b) the scalar particle. 

of the spin. The obtained result is quite expectable if we recall that in the case of a weak 
field, the 5-well itself can bind a charged particle only if it is in the s state. 

The dependencies of ln(| J|/M) on the transverse coordinate in the planes with different 
z for the strong magnetic field h = 100 demonstrating the spatial distribution of the 
probability currents are shown in Fig. 9a for the electron with the spin directed opposite 
the field and in Fig. 9b for the scalar particle. Exact formula (24) is used in the calculation. 

Figure 10 shows the comparison of the behavior of ln(| J|/M) for the scalar particle 
with that for the electron with its different spin orientations in the case of the strong field 
and the plane z = 0. In particular, it can be seen from Fig. 8 that the orbital and spin 
currents in the strong field become comparable in the vicinity of the 5-center. As should 
be expected, in the case of the electron with the spin oriented opposite the magnetic field, 
the current decreases more rapidly as the distance from the center increases, i.e., this 
particle motion is localized. The range where the orbital and spin currents turn out to be 
comparable is probably most interesting for studies. 
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Figure 10: The dependence of ln(|J|/M) on the transverse coordinate obtained using 
exact formula (24) in the case of the plane z = and the strong field h = 100 for different 
spin orientations (along the field and opposite the field) and the zero spin. 
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